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Abstract 

The equilibrium of a cylindrical plasma with purely poloidal mass flow and 
cross section of arbitrary shape is investigated within the framework of the 
ideal MHD theory. For the system under consideration it is shown that only 
incompressible flows are possible and, conscequently, the general two dimen- 
sional flow equilibrium equations reduce to a single second-order quasilinear 
partial differential equation for the poloidal magnetic flux function if), in which 
four profile functionals of ip appear. Apart from a singularity occuring when 
the modulus of Mach number associated with the Alfven velocity for the 
poloidal magnetic field is unity, this equation is always elliptic and permits 
the construction of several classes of analytic solutions. Specific exact equlib- 
ria for a plasma confined within a perfectly conducting circular cylindrical 
boundary and having i) a flat current density and ii) a peaked current density 



-mail: gthroum@cc.uoi.gr 
-mail: gpantis@cc.uoi.gr 

be published in Plasma Physics and Controlled Fusion 



1 



are obtained and studied. 
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I. INTRODUCTION 



Plasma rotation of the order of sound speed has been observed in many tokamaks heated 
by neutral beams, such as JET [[]J and TFTR [0. Early reported data relates to purely 
toroidal rotation 0; more recent measurements show that poloidal rotation can occur as 
well 0]. It has also been found that poloidal flows in the outer region of tokamaks driven 
by strong electric fields play an important role in the L-H transition (e.g. |5|,[]]). 

Although theoretical studies of flowing plasmas began in the mid fifties J7|- [[|, since the 
early eighties there has been an increasing interest (motivated from the above mentioned 
observations) mainly in the investigation of equilibrium and stability of plasmas with mass 
flow |10|]- f22j. In particular, the symmetric equilibrium states of a plasma with mass flow 
in a two dimensional geometry are governed by a second-order partial semilinear differential 
equation for the poloidal magnetic flux function -0, which contains five free surface quanti- 
ties (i.e. quantities solely dependent on ip) in conjuction with a nonlinear algebraic Bernouli 
equation (e.g. |10|- fl7|l ). Unlike in static equilibria, the above mentioned differential equa- 
tion is not always elliptic; there are three critical values of the modulus of Mach number 
\M\ of the poloidal flow with respect to the poloidal magnetic field, at which the type of 
this equation changes, i.e. it becomes alternatively elliptic and hyperbolic. To solve the 
equilibrium problem in the two elliptic regions, several computer codes have been developed 
(e.g. |l5|-|r7|]). For purely toroidal flow the partial differential equation becomes elliptic and 
can be solved analytically; analytic solutions were obtained in Refs. |23| - [ 25(1 - Owing to 
the complexity of the equilibrium equations, the construction of analytic solutions in the 
case of arbitrary flows (toroidal and poloidal) is difficult. To our best knowledge the only 



analytic solution of this kind is that obtained by Mashke and Perrin |26| in the specific case 
of small ratio of poloidal to toroidal magnetic field and small (3 defined as the ratio of the 
kinetic energy to the magnetic energy. 

In the present paper we study the equilibrium of a cylindrical plasma with purely poloidal 
mass flow in the framework of ideal MHD theory. Several physics aspects of a similar system 
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has been often considered in the literature, e.g. i) the numerical study of Ref. concerning 



the current sheet formation in a cylindrical plasma with rectangular cross section and purely 
poloidal flow and ii) the analytic study of Ref. on nonlinear interactions of tearing modes 
in plane geometry in the presence of sheared poloidal flow. The plasma cross section of the 
system under consideration is arbitrary and therefore the configuration is two-dimensional. 
No restriction is imposed either on the components of the magnetic field or j3. We show 
that only incompressible equilibrium flows are possible and that several classes of analytic 
equilibrium solutions can be obtained. 

The equilibrium equations of the system under consideration are presented in Section 
II. In Section III we derive exact equilibrium solutions. Our conclusions are summarized in 
Section IV. 



II. EQUILIBRIUM EQUATIONS 



The system under consideration is a cylindrical plasma with mass flow and arbitrary 
cross section shape. For this configuration convenient coordinates are £, rj and z with unit 
basis vectors e^, e^, e z , where e z is parallel to the axis of symmetry and £, i] are generalized 
coordinates pertaining to the poloidal cross section. Specific examples are circular and 
elliptical cylindrical coordinates. The equilibrium quantities do not depend on z. The 
equations governing the equilibrium states can be obtained from the general two dimensional 



equilibrium equations with flow, e.g. Refs. [11|, [[12], fL5 |. In particular, the components of 
the momentum conservation equation 



p(v • V)v = j x B - VP 



along the magnetic field and perpendicular to a magnetic surface, respectively, are put in 
the form 



B 



vl \VF\ 2 \ VP 



V — + 
I 2 2p 2 



+ 



P 







(1) 



and 



V 



P |V^| 2 



+ 



F"F'|V^| 2 B Z VB Z -V^ 



VP" 

p 



0. 



(2) 



Here, ipi^tf), F(£,rj), B z (£,rj) and v z (£,rj) are stream functions in terms of wich the diver- 
gence free fields, i.e. the magnetic field B, the current density j and the mass flow pv, are 
expressed as 



B = B z e z + e z x V^, 



(3) 



j = V 2 ipe z -e z x VB Z 



(4) 



and 



pv = pv z e z + e z x VF; 



(5) 



the function F and the electrostatic potential $ are surface quantities satisfying the equa- 
tions 



B z - F'(VH = XW>) 



(6) 



(X also is a surface quantity.) and 
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z F'{iP)-v z = &{iP); 



P 



(7) 



the prime denotes derivative with respect to i/j; SI units are used and the vacum magnetic 
permeability /z is set to unity. It is noted here that a particular equation of state, e.g. 
either isentropic or isothermal magnetic surfaces, has not been employed in order to express 
in Eq. (fj) the pressure gradient P' in terms of surface quantities; thus, Eq. (|2]) holds for 
any equation of state. 

We now restrict attention to purely poloidal flow and therefore set v z = 0. Then, Eq. 
(§) implies that the function B z is a surface quantity, i.e. B z = B z (ip). From equation (|7|) 
it follows then that 



P = pW, (8) 

and therefore the mass density and the magnetic field share also the same flux surfaces. 
Using Eq. (§), the mass conservation equation V ■ (pv) = yields V ■ v = 0. Thus, it 
turns out that for a cylindrical plasma and purely poloidal mass flow only incompressible 
equilibrium flows are possible. 

With the use of Eq. @, Eq. (]!]) can now be readily integrated, yielding an expression 
for the pressure, i.e. 

p = m)-—m 2 . (9) 

We note here that, unlike in static equilibria, in the presence of flow magnetic surfaces do 
not coincide with isobaric surfaces because the momentum balance equation implies that 
B ■ VP in general differs from zero. In this respect, the term P s (ip) is the static part of 
the pressure which does not vanish when F' is set to zero; thus, the second term makes the 
isobaric surfaces to detach from magnetic surfaces. 

Before proceeding to obtain the final form of the equilibrium equation we examine a 
particular case, i.e. when the functions F and p satisfy the relation 

^ = 1. (10) 

Then, the second-order differential operator in Eq. (Q) is singular and since F and p are 
surface quantities this singularity occurs uniformly over a magnetic surface. On the basis 
of Eq. ([|) for pv and the definitions v 2 Ap = - — for the Alfven velocity associated with 

2 

the poloidal magnetic field and the Mach number M 2 = , Eq. (pi]) can be written as 

v Ap 

( F') 2 

M 2 = 1. Assuming now - ^ ^ 1, inserting Eq. (|S|) into Eq. (^) and setting v z = 0, Eq. 
(0) after some straightforward algebra, reduces to the elliptic differential equation 
{F'f 



1 - 



W + 7 (y^ - F ") I v ^l 2 + ( p ° + t ) = o- ( n ) 



p 

The absence of any hyperbolic regime in Eq. (|TTJ) can be understood by noting that, as is well 
known from the gass dynamics, the flow must be compressible to allow the governing partial 
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differential equation to depart from ellipticity. It may also be noted that in axisymmetric 
equilibria with arbitrary flows, e.g. Eq. (0), \M\ = 1 is not a transition point. Finally, we 
note that replacing the stream functions by the corresponding field quantities, Eq. ( |TTD is 
put in the form 

d ( |B| 2 \ B 2 v 2 , . 

crr[ p+ -i) + T -'7 = °' (12) 

the terms of which have a readily unterstandable physical meaning. 

III. EXACT SOLUTIONS 

In this section, Eq. fllTl) is solved for a specific choice of the flux functions p{i]S) and 
F(ip), i.e. when these functions satisfy the relation 

p' F" , . 

7 = V (13) 



This ansatz makes the second term in Eq. (JTTj) to vanish. Eq. (|T3| ) is an ordinary differential 
equation for if), which can easily be integrated. Its solution is 

= M 2 , (14) 

where M 2 =const. ^ 1. The choice ( ]TB| ) leads therefore to flows with constant Mach 
numbers. Solution (|i~4"D makes the factor multiplying the second-order differential operator 
in Eq. (fTID constant; thus, Eq. ([0]) reduces to 

v ^ + T^( F « + f)'-°- < 15 > 

This is similar in form to the equation governing static equilibria; the only explicit reminis- 



cent of flow is the presence of M c . Eq. (|15D, which is still a semilinear partial differential 

B 2 B 2 
equation, can be linearized for several choices of P s H — e -§- a ) f° r P s H — linear in 

x/j becomes an inhomogeneous Laplace equation with constant inhomogeneouity and b) for 
B 2 

P s H — tj- quadratic in ip becomes i) a homogeneous Helmholtz equation for M 2 > 1 and 
ii) a homogeneous diffusion equation for < M 2 < 1. Several classes of analytic solutions 



can therefore be derived. In the following, specific exact equilibria will be constructed for a 
circular cylindrical plasma having i) a flat current density profile and ii) a peaked current 
density profile. 



A. Flat Current Density Profile 



To describe a plasma which is confined within a circular infinitely conducting cylinder 

of radius r , the coordinates £, rj, z are specified to be the circular cylindrical coordinates 

r, 9, z with unit basis vectors e r , e e , e z . The equilibrium quantities depend just on r. To 

( B 2 \' 

obtain equilibria having a flat current density profile, we employ the ansatz ( P s + -if- J = 
const. The simplest non trivial solution of Eq. (|T5| ) is written then as 



ip oc r 



(16) 



where r = -f^. Solution ( |16D yields a singly peaked parabolic pressure profile 



P(p) = F(0)(l-r 2 ), 



(17) 



a flat "toroidal" current density 
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and 
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F' 
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(22) 



Here, a is a parameter which, together with M c , describes the magnetic properties of 
the plasma, i.e. the plasma is diamagnetic for ^- — ° L ^ , i ^ < 1 an d paramagnetic for 

(l — ~\<p ) > ^' van ^ snm S now {Mc — 0) the equilibrium reduces to the exact ax- 
isymmetric equlibrium solution, in the limit of infinite aspect ratio, which was suggested by 
Shafranov and has been extensively used by Solov'ev [^9j. In Eqs. (^) and (p2|), the 
density p remains a free function of r and therefore a variety of velocity and electric field 
profiles can be derived. As an example, we consider an equilibrium having a singly peaked 
density profile: 

-- p{0 \k- P r, (23) 



k" 



where v > and k > 1 . The parameter k here is necessary in order to make p non- vanishing 
on plasma surface (p = 1). (Note that p appears in the denominator of Eqs. (|21|) and (p2[)). 
This agrees qualitatively with experimental results |6,27], according to which finite densities 



and very low temperatures were measured in the edge region of plasmas with mass flow. 
For the density profile (E3p, the velocity profile becomes hollow and the temperature profile, 
consistent with Eq. flTT|), is given by 

r = T(0)**£^. (24) 

B. Peaked current density profile 

For simplicity we ignove here diamagnetic effects, viz. we consider (3 P = 1 equilibria. 
(The construction of f3 p ^ 1 equilibria is, however, possible). Accordingly, for B z = B = 
const, and 

F s a(l-M c 2 )^ 2 



the simplest solution of Eq. ([iq) is 

^ = #))J (C). (25) 
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Jo and J\ below are zeroth- and first-order Bessel functions, respectively; £ = , where 
A is the first zero of J . Solution (ESI) yields a peaked pressure profile 



P = P.(0) [(1 - M c 2 ) J (0 + M 2 (J 2 (A ) - J 2 (C))J , (26) 
a singly peaked current density profile 

J z = JMUQ, (27) 



Be = -^mMO, (28) 



±|MjM!)« e , (29) 

^0 v 7 ^ 



and 



E^IMjMW*,, (30) 

Here, P s (0) = P(0)|m c =o is the static pressure at the plasma center. The above equilibrium 
has the following characteristics: 

1. Both the pressure and the current density profiles vanish on the plasma surface. 

2. The magnetic field is independent of flow. 

We note that the requirement P(0) > imposes a restriction on the values of Mach number, 
i.e. M 2 < [1 — Ji(Aq) 2 ] . By inspection of Eq. (|26| ) it turns out that as flow increases, 



i.e. when \M C \ takes higher values, the pressure P takes lower values and its profile becomes 
flatter; the negative term takes therefore higher values. On the other hand, as |v| 



2 

increases, the flow term —py- (for a given density profile) becomes more negative. Thus, 
HP ?; 2 

the quantity ^f- — pjr becomes flow independent and therefore guarantees the equilibrium 
equation (0), in which the magnetic field terms (for the solution (|25[) ) are independent of 
flow. 
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IV. CONCLUSIONS 



Using the general two dimensional ideal MHD equlibrium equations with flow, we inves- 
tigated the stationary states of a cylindrical plasma with purely poloidal flow and arbitrary 
cross-section shape and showed that the flows are necessarily incompressible. Consequently: 
(i) The pressure profile is expressed in terms of a static pressure surface quantity P s and 
a flow dependent term not remaining constant on a magnetic surface (Eq. (0)), which is 
consistent with the fact that for equilibria with mass flow the magnetic surfaces are detached 
from the isobaric surfaces, (ii) The flow equations reduce to a single second-order quasilinear 
partial differential equation for ip, (Eq. fll!] )), which contains four surface quantities, i.e. 
P 8 (if > ), the flow function F(if)), the "toroidal" magnetic field B z {ijj) and the density p(i]j). 
One has to specify these quantities together with the boundary conditions for i/j in order 
to determine the equilibrium completely. A singularity occurs in the differential equation 
when the modulus of the Mach number \M\, associated with the Alfven velocity for the 
poloidal magnetic field, is unity. Apart from this case, this equation is always elliptic. This 
is consistent with the fact that the existence of hyperbolic regimes in the governing equilib- 
rium differential equation is possible only in compressible fluids. For velocities satisfying the 
relation \M\ — const. 7^ 1, the equilibrium equation takes a simpler form (Eq. (0)), which 
is amendable to several classes of analytic solutions. We constructed and studied particu- 
lar exact equilibria for a plasma confined within a circular cylindrical perfectly conducting 
boundary having i) a flat current density profile and ii) a peaked current density. 

The analytic equilibria with poloidal flow we constucted in the present paper can be 
the basis of stability and transport investigations, which would be of relevance to magnetic 
confinement systems. In particular, they may help in understanding the physics involved in 
the transition from the low confinement mode to the high confinement mode in tokamaks, 
which is crucially related to poloidal flow. The derivation of new analytic solutions in more 
realistic situations, e.g. axisymmetric equilibria, hellicaly symmetric equilibria (even in 
specific cases such as in rotating plasmas with a single velocity component or/and small /?), 



11 



might help to easily get physical insight of plasmas with flows. Solutions of this kind might 
also be used for checking relevant equilibrium codes. 
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